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From:  A.  L.  Treaster,  P.  P.  Jacobs,  Jr.*  and  G.  B.  Gurney 

Subject:  Sidewall  Boundary  Layer  Corrections  in  Subsonic,  Two-Dimensional 

Airfoil/Hydrofoil  Testing** 

Abstract:  Historically,  two-dimensional  airfoil  or  hydrofoil  section 

characteristics  have  been  obtained  by  measuring  individually  the  lift,  drag 
and  pitching  moment  by  the  most  accurate  technique  available.  The  use  of 
force  balances  to  measure  the  three  quantities  simultaneously  has  met  with 
only  partial  success.  Although  the  lift  and  pitching  moment  data  have  usually 
been  acceptable,  the  drag  data  have  varied  by  as  much  as  an  order  of  magnitude 
from  previous  reference  data.  To  investigate  the  parameters  which  influence 
two-dimensional  force  measurements,  an  experimental  program  was  conducted 
in  the  subsonic  wind  tunnel  of  the  Applied  Research  Laboratory  at 
The  Pennsylvania  State  University.  From  the  results  of  this  test  program, 
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the  sidewall  boundary  layer  was  identified  as  the  primary  factor  contributing 
to  the  erroneous  drag  measurements.  A  correction  procedure  which  is  based  on 
the  airfoil/hydrofoil  geometry,  the  flow  environment  and  the  measured  data  was 
developed.  Corrected  data  from  the  subject  test  program  and  from  similar 
programs  in  other  experimental  facilities  for  both  symmetrical  and  cambered 
sections  are  in  good  agreement  with  the  reference  data  in  all  cases. 
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Nomenclature 

2 

aspect  ratio  =  s  /sc 

exponent  in  Eq .  (11) 

airfoil  chord  length 

sectional  drag  coefficient  =  D/Cq^s) 

sectional  drag  coefficient  at  c^  =  0.0 

sectional  lift  coefficient  =  L/(q  s) 

oo 

local  slope  of  the  c^  vs  a  curve 

slope  of  the  linear  portion  of  the  c^  vs  a  curve 
(usually  evaluated  in  the  c^  =  0.0  region) 

the  drag  force 

Hawthorne's  approximation  of  the  energy  in  a  secondary  flow 
[Eq.  (1)] 

defined  by  Eq .  (2)  and  Figure  10 
functional  operators  used  in  developing  Eq.  (12) 
proportionality  constants  used  in  developing  Eq .  (12) 
the  lift  force 

defined  by  Eq .  (3) 
atmospheric  pressure 
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Pg  =  static  pressure 

PTOTAL  =  total  pressure 

2 

q  =  local  dynamic  pressure  =  l/2pu 

2 

q  =  free  stream  dynamic  pressure  =  l/2pV 

00  CD 

Re  =  Reynolds  number  =  V^c/v 

s  =  airfoil  span 

t  =  maximum  airfoil  thickness 

u  =  local  velocity 

=  free  stream  velocity 

x  =  distance  parallel  to  test  section  centerline  measured  from 

the  leading  edge  of  the  two-dimensional  test  chamber 

a  =  airfoil  angle  of  attack 

Ac,  =  required  correction  to  c, 
a  a 

6  =  test  section  wall  boundary  layer  thickness  at  the  centerline 

of  the  balance  shaft 

* 

6  =  test  section  wall  boundary  layer  displacement  thickness  at 

the  centerline  of  the  shaft 

p  =  mass  density  of  the  fluid 


v 


kinematic  viscosity  of  the  fluid 
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Introduction 

The  use  of  a  force  balance  to  measure  simultaneously  the  sectional 
characteristics  of  an  airfoil  or  hydrofoil  that  spans  a  rectangular  test 
section  has  met  with  only  partial  success.  Although  the  lift  and  pitching 
moment  data  have  usually  been  acceptable,  the  drag  data  (with  the  traditional 
corrections  applied)  have  varied  by  as  much  as  an  order  of  magnitude  from 
established  reference  data.  Based  on  studies  conducted  in  the  subsonic  wind 
tunnel  of  the  Applied  Research  Laboratory  at  The  Pennsylvania  State  University 
(ARL/PSU)  the  sidewall  boundary  layer  was  identified  as  the  primary  factor 
contributing  to  these  erroneous  drag  measurements.  Presented  herein  is  an 
empirically  derived  method  to  correct  for  the  presence  of  the  sidewall 
boundary  layer  in  subsonic  two-dimensional  airfoil/hydrofoil  testing  with  a 
mechanical  force  balance. 


Background 

A  review  of  the  literature  indicated  that  most  of  the  currently  used  NACA 
airfoil  section  data  were  measured  in  wind  tunnels  during  the  1930's  and 
1940' s  at  both  the  Langley  and  Ames  Research  Centers.  The  majority  of  these 
data  were  acquired  by  measuring  individually  the  lift,  drag  and  pitching 
moment  by  the  most  accurate  means  available.  Generally,  this  required  either 
the  measurement  of  the  pressure  distributions  on  the  ceiling  and  floor  of  the 
test  section  or  the  use  of  a  force  balance  to  obtain  lift.  Drag  data  were 
obtained  from  either  wake  surveys  or  via  surface  pressure  distributions;  a 
torsional  balance  was  usually  used  to  measure  the  pitching  moment.  Typical 
of  these  measurement  programs  were  those  conducted  by  Loftln  and  Smith! 
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In  wind  tunnel  testing,  the  sidewall  boundary  layer  can  sometimes  be 
managed  by  blowing  or  suction  techniques.  More  recently  computational  * 

procedures,  such  as  those  by  Barnwell2,  Sewall^,  and  Kemp  and  Adcock4,  have 
been  applied  to  the  problem  for  selected  flow  fields. 

As  would  be  expected,  the  sidewall  boundary  layer  problem  is  also  present  * 

in  two-dimensional  hydrofoil  testing.  Kermeen’’  and  Daily^  both  used 
mechanical  force  balances  at  The  California  Institute  of  Technology  (CIT) 
to  measure  forces  on  hydrofoils.  In  general,  their  lift  and  pitching  * 

moment  data  were  in  good  agreement  with  previous  measurements.  Although 
their  drag  measurements  agreed  well  with  the  available  reference  data  in  the 
low  angle  of  attack  range,  at  larger  values  of  a  their  measurements  were  high  * 

by  as  much  as  a  factor  of  two.  In  the  early  1970's,  researchers  at  ARL/PSU, 
in  an  effort  to  develop  new  propulsor  blade  design  criteria,  used  a  three- 
component,  mechanical  force  balance  in  an  attempt  to  measure  the  two- 
dimensional  sectional  characteristics  of  a  hydrofoil  that  spanned  the 
rectangular  test  section  of  the  ARL/PSU  12  in.  (304.8  mm)  cavitation  tunnel2. 

► 

With  the  exception  of  the  pitching  moment  characteristics,  the  data  were  in 
disagreement  with  available  reference  data  (later  reorientation  of  the 
sensing  elements  solved  the  lift  problem). 

» 

Although  air  and  water  are  both  fluids,  fundamental  differences  exist  in 
the  testing  requirements  between  the  two  media.  Testing  airfoil  shapes  in 
water  (hydrofoils)  introduces  additional  problems  such  as  the  handling  of 

► 

larger  gross  forces  and  waterproofing  requirements.  However,  the  major 
concern  in  water  tunnel  measurements  is  the  cavitation  phenomenon.  A  hydro¬ 
foil  may  operate  in  any  or  all  of  three  different  flow  regimes;  namely, 

> 

ful Ly  wetted  flow,  partially  cavitating  flow  or  fully  cavitating  flow. 
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Testing  in  fully  wetted  flow  differs  little  from  low-speed  wind  tunnel 
testing,  and  it  is  this  regime  that  was  of  primary  concern  in  the  original 
ARL/PSU  water  tunnel  test  program. 

In  a  water  tunnel,  placing  pressure  taps  on  a  hydrofoil  surface  or  on 
the  test  section  walls  can  cause  premature  cavitation  and  result  in  erroneous 
pressure  measurements.  At  certain  flow  conditions,  pressure  probes  in  the 
hydrofoil  wake  are  also  subject  to  cavitation  problems.  For  these  reasons 
water  tunnel  force  measurements  are  best  performed  by  a  mechanical  force 
balance.  With  such  a  balance,  forces  can  be  measured  directly  without  marring 
the  model's  surface.  As  previously  discussed,  mechanical  balances  are  not 
totally  free  of  problems.  Balances  measure  all  forces  applied  to  a  model; 
and  if  forces  occur  on  a  model  which  are  not  those  associated  with  two- 
dimensional  flow,  the  balance  will  also  measure  them.  A  correction  procedure 
is  then  required  to  reduce  this  measured  total  force  to  a  two-dimensional 
force.  However,  the  traditional  corrections  discussed  by  Pope^  and  Allen 
and  Vincenti^;  namely,  solid  blockage,  wake  blockage,  lift  effect  and 
horizontal  buoyancy  are  not  sufficient  to  satisfactorily  correct  the  drag 
data.  Thus,  the  diagnostic  test  program  to  identify  additional  correction 
procedures  was  conducted  in  the  ARL/PSU  subsonic  wind  tunnel^.  This  facility 
provided  an  environment  in  which  the  aerodynamic  and  geometric  parameters 
common  to  both  airfoil  and  hydrofoil  testing  could  be  easily  and  economically 
varied.  The  results  of  this  investigation  conducted  by  Jacobs^  verified 
that  the  two-dimensionality  of  the  flow  was  being  contaminated  by  the  inter¬ 
action  of  the  airfoil  with  the  sidewall  boundary  layer  and  permitted  the 
formulation  by  Jacobs  of  an  empirical  correction  procedure  which  is  summarized 
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here  In  Eq.  (12).  When  this  correction  procedure  Is  applied  to  the  ARL/PSU- 
measured  drag  data  and  to  additional  drag  data  measured  by  Ward^  at  C1T, 
the  results  are  in  good  agreement  with  the  reference  data. 

The  specific  details  of  the  entire  wind  tunnel  test  program  are 
documented  in  the  report  by  Jacobs^.  Discussed  in  the  remainder  of  this 
paper  are  the  portions  of  the  test  program  relevant  to  the  lift  and  drag 
measurements,  the  development  of  the  correction  procedure  and  the 
application  of  the  correction  procedure  to  existing  experimental  data. 

ARL/PSU  Wind  Tunnel  Test  Program 

Test  Facility  and  Experimental  Hardware 

The  ARL/PSU  subsonic  wind  tunnel  is  a  closed  circuit,  closed  jet  wind 
tunnel  with  an  octagonal  test  section  which  is  4.0  ft  (1.219  m)  across  the 
flats  and  is  16.0  ft  (4.877  m)  long.  The  test  section  velocity  can  be  varied 
continuously  up  to  120.0  fps  (36.576  m/sec).  Honeycomb  and  screens  used  in 
the  settling  section  reduce  the  turbulence  level  in  the  test  section  to  less 
than  0.10  percent  of  the  free-stream  velocity  at  80.0  fps  (24.384  m/sec). 

For  this  test  program,  two  4,0  ft  x  8.0  ft  (1.219  m  x  2.438  m)  wooden  panels 
were  mounted  vertically  18.375  in.  (466.725  mm)  apart  to  create  a  two- 
dimensional  test  section  as  shown  in  the  balance  installation  drawing, 

Fig.  I. 

Two  NACA  0012  airfoils  of  aspect  ratios  1.02  (c  =  18.0  in.  (457.2  mm), 
and  s  =  18.375  in.  (466.7  mm))  and  2.04  fc  =9.0  in.  (228.6  mm);  and 
s  =  18.375  in.  (466.7  mm)]  were  fabricated.  When  installed  in  the  test 
chamber,  the  airfoil  was  attached  at  its  midchord  to  the  balance  by  a 
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B 

■ 

m 

a 


spanwlse  Internal  shaft.  The  balance  shaft  was  located  midway  between  the 
floor  and  ceiling  of  the  test  section  and  28.0  in.  (711.2  mm)  downstream  from 
the  leading  edge  of  the  wooden  panels. 

Instrumentation 

Because  measurement  of  the  pitching  moment  had  not  been  a  problem  in  the 
past,  lift  and  drag  forces  only  were  measured  by  a  two-component  balance 
that  is  sketched  in  Fig.  1.  This  balance  used  the  compact  strain-gaged 

1  O 

tension-member  concept  developed  by  Gurney1  .  The  balance  rotated  with 
the  airfoil  and  sensed  forces  normal  to  and  parallel  with  the  chordline. 

For  the  force  measurements  the  reference  velocity,  ,  was  measured  by 
a  0.25  in.  (6.35  mm)  diameter  pitot-static  probe  that  was  located  on  the 
tunnel  floor  midway  between  the  sidewalls  and  in  the  same  vertical  plane  as 
the  midchord  of  the  airfoil.  The  probe  tip  was  above  the  floor  boundary 
1  aye  r . 

Wake  traverses  to  evaluate  the  sectional  drag  coefficient  were  conducted 
with  a  0.125  in.  (3.175  mm)  diameter  kiel  probe  that  was  located  midway 
between  the  sidewalls  and  in  a  plane  one  chord  length  downstream  from  the 
trailing  edge  of  the  airfoil.  For  these  tests,  the  reference  pitot-static 
probe  was  located  in  the  traverse  plane  midway  between  the  floor  and  the 
centerline  of  the  test  section,  Fig.  2.  The  same  kiel  probe  was  used  to 
make  sidewall  boundary  layer  measurements  at  the  balance  shaft  location.  For 
these  measurements,  the  reference  pitot-static  probe  in  the  tunnel  floor 


was  used 
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To  obtain  horizontal  buoyancy  corrections,  the  sidewall  static  pressure 
gradient  was  measured  by  four  static  pressure  taps  along  the  horizontal 
centerline  of  the  sidewall.  The  angle  of  attack  was  measured  with  a 
gunners  quadrant  in  conjunction  with  an  accurately  machined  airfoil  template. 

Measurements  and  Results 

Test  Section  Characteristics 

The  flow  characteristics  of  the  two-dimensional  test  section  were 
established  by  a  series  of  preliminary  tests  prior  to  the  installation  of  the 
airfoil.  Flow  uniformity  was  verified  for  the  region  outside  of  the  influence 
of  the  four  test  section  boundary  layers  by  kiel  probe  surveys.  The 
longitudinal  static  pressure  gradient  [  (dP  /dx)/l/2pV^J  was  measured  to  be 
0.01236  ft  *  (0.0406  m  1).  The  sidewall  boundary  layer  at  the  balance  shaft 
location  was  measured  at  several  Reynolds  numbers;  the  resulting  data  are 
presented  in  Fig.  3. 

Establishing  a  Reference  Data  Base 

Because  no  reputable  drag  polars  could  be  located  at  the  target  Reynolds 
number  (Re  =  330,000),  new  baseline  data  were  measured  as  a  part  of  the  wind 
tunnel  test  program.  The  proven  NACA  approach  of  using  separate  lift  and  drag 
measurements  with  modern  instrumentation  was  utilized  with  the  9.0  in.  (228.6) 
chord  airfoil.  The  lift  data  were  measured  by  the  force  balance,  whereas  the 
drag  data  were  obtained  from  momentum  principles  applied  to  downstream  wake 
traverses  at  the  midspan  of  the  airfoil.  In  Fig.  4  the  resulting  lift  data 
are  shown  in  comparisons  with  the  Loftin  and  Smith^  measurements  at  a  Re  of 
700,000  and  with  computational  data  from  Ohio  State  University^  (0SU)  at 
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Re  =  330,000.  For  completeness  the  balance-measured  data  with  18.0  in. 

(457.2  mm)  chord  airfoil  are  also  shown.  These  data  were  judged  to  be  in 
satisfactory  agreement. 

The  resulting  drag  data  are  shown  by  the  solid  "diamonds"  in  Fig.  5.  To 
establish  a  trend  with  Re,  data  from  other  sources  are  also  included  in 
Fig.  5.  The  current  data  fit  in  well  with  the  observed  Re  variation.  It 
should  be  noted  that  the  only  reference  data  at  the  target  Re  are  the 


computational  data  from  OSU^  and  the  high  turbulence  level  data  of  Jacobs 
and  Sherman^.  The  OSU  computational  data  predicts  the  onset  of  separation 
but  does  not  include  the  associated  effects  on  c^  and  in  the  computations 
and,  thus,  underpredict  c^  at  the  higher  c^  values.  For  these  reasons,  the 
ARL/PSU  sectional  characteristics  measured  by  the  combined  balance-wake 
traverse  approach  were  used  as  the  "reference  data”  for  the  remainder  of  the 
study . 


Balance  Measurements 

The  results  of  the  balance  measurements  are  shown  in  Figs.  4  and  6. 

In  Fig.  6,  the  discrepancy,  Ac^  ,  between  the  balance-measured  drag  and 
the  reference  data  is  evident. 

What  was  the  source  of  error  in  the  balance-measured  drag  data?  Three 
possible  sources  were  considered:  (1)  end  gap  effects  between  the  airfoil 
tip  and  the  adjacent  sidewall,  (2)  drag  on  the  portion  of  the  balance  shaft 
between  the  sidewall  and  the  model,  and  (3)  contamination  of  the  two- 
dimensional  flow  by  an  interaction  of  the  airfoil  and  the  sidewall  boundary 
layer.  The  effect  of  the  end  gap  between  the  airfoil  tip  and  the  test 
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chamber  wall  was  first  investigated.  Lift  and  drag  data  were  measured  at 
a  constant  angle  of  attack  while  varying  the  end  gap  from  0.001  in.  to 
0.010  in.  (0.025  ram  to  0.254  mm).  For  this  range  of  end  gap,  no  significant 
change  in  c^  was  measured  and  less  than  a  1 .0  percent  change  in  c ^  was 
recorded.  Parkin  and  Kermeen^  reported  similar  results,  i.e.,  if  the  end 
gap  is  sufficiently  small,  viscous  forces  predominate  and  the  effect  of  the 
end  gap  is  negligible. 

To  evaluate  the  drag  force  on  the  portion  of  the  balance  shaft  exposed 
to  the  flow  a  stub  spindle  was  fabricated.  The  stub  spindle  was  mounted  in 
the  balance  and  extended  0.002  in.  (0.508  mm)  into  the  flow.  This  was  the 
typical  operating  clearance  at  ‘.he  balance  end  of  the  airfoil.  The  9.0  in. 
(228  .5  mm)  airfoil  was  mounted  on  the  opposite  wall  and  was  maintained 
at  a  minimum  distance  from  the  spindle.  Under  these  conditions,  the 
flow  in  the  vicinity  of  the  model-sidewall  intersection  was  closely 
duplicated  and  the  balance  measured  only  the  forces  on  the  stub  spindle. 

The  effect  was  negligible. 

Thus,  only  the  contamination  of  the  two-dimensional  flow  by  the  inter¬ 
action  of  the  airfoil  with  the  sidewall  boundary  layer  remained  as  a 
postulated  cause  of  the  erroneous  drag  measurement.  Evidence  of  such  an 
interaction  was  observed  during  the  wake  measurement  phase  where  a  secondary 
wake  was  measured  near  the  sidewall.  Fig.  7.  As  previously  discussed,  the 
removal  of  the  sidewall  boundary  layer  is  not  practical  in  water  tunnel 
applications  so  that  the  alternate  approach  of  developing  a  correction 
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procedure  was  chosen. 
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Development  of  the  Correction  Procedure 

Shown  in  Fig.  6  is  the  variation  between  balance-measured  drag  data  and 

the  ARL/PSU  reference  two-dimensional  section  characteristics  for  the 

NACA  0012  airfoil.  Because  this  difference  is  typical  only  of  balance- 

measured  data  in  which  the  entire  force  exerted  or  the  airfoil/hydrofoil  is 

measured,  it  was  assumed  that  this  Increment  in  drag,  Ac,  ,  was  due  to 

d 

three-dimensional  flow  effects  on  the  model.  Such  three-dimensional  flow 
effects  can  be  generated  when  a  strut  intersects  a  flat  surface  in  the 
presence  of  a  nonuniform  flow.  The  resulting  secondary  flow  —  the  so-called 
horseshoe  vortex.  Fig.  8  —  engulfs  the  strut-wall  intersection  and  produces 
a  region  of  contaminated  two-dimensional  flow.  This  type  of  secondary  flow 
can  be  generated  in  airfoil/hydrofoil  testing  when  an  airfoil  or  hydrofoil 
that  spans  the  test  section  intersects  the  test  section  wall  in  the  presence 
of  the  sidewall  boundary  layer. 

A  study  of  this  problem  using  flow  visualization  techniques  was  recently 
completed  by  Barber^  in  which  he  investigated  the  additional  drag  that  is 
created  by  a  strut  protruding  from  a  wall  as  a  function  of  the  incoming 
boundary  layer  thicknesses.  He  found  that  the  size  of  the  horseshoe  vortex 
varied  directly  with  the  thickness  of  the  incoming  boundary  layer.  He  also 
found  that  the  portion  of  the  airfoil  where  flow  separation  occurred  varied 
inversely  with  the  size  of  the  horseshoe  vortex.  He  concluded  that  with  a 
large  horseshoe  vortex,  viscous  effects  caused  high  energy  fluid  to  be 
entrained  in  the  corner  where  the  airfoil  trailing  edge  and  wall  intersects 
as  shown  in  ”i  ».  q .  This  influx  of  high  energy  fluid  enables  the  flow  to 
withstan  '  h u  ter  the  adverse  pressure  gradient  existing  in  the  corner  and. 
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consequently,  retards  flow  separation.  ,4s  illustrated  in  Fig.  9,  a  thin 
vortex  is  not  able  to  entrain  as  much  of  the  high  energy  fluid  and  a  larger 
separated  zone  exists. 

Hawthorne^  derived  the  following  expression  for  the  energy  in  secondary 
* 

flows,  D  ,  created  by  strut-wall  intersections: 
e 

144V^c“(t/c)4f (n) 

D  =  - =~  ,  O) 

25  [  1  +  <l/2)(t/c)Z] 

where 


f  (n)  = 


1  +  n 


1  {iL_ 

2  it  4n 


n  +  1 


1  -  n 


(1  +  nZ) 


T  lo*e  n  + 


l  +  n 


2>  “ft 


(2) 


and 


n  =  4  [  1  +  (l/2)  (t/c  )  ]/(15tt6  /c)  .  (3) 

Hawthorne's  relationship  between  f(n)  and  boundary  layer  displacement 

thickness  (6  /c)  is  shown  in  Fig.  10.  These  data  are  for  a  bicusped  strut 

profile  in  an  exponential  boundary  layer  with  strut  thickness-to-chord  ratios 

* 

of  .05  and  .25.  Hawthorne's  figures  show  that  f(n)  increases  with  6  /c  to 
a  maximum  value  at  6  / c  =0.1  Equation  (l)  states  that  the  energy  in  these 
secondary  flows  is  proportional  to  airfoil  thickness  to  the  fourth  power  and 

'k 

reaches  a  maximum  when  6  /c  is  approximately  0.1.  Although  the  theory  does 
not  hold  for  all  airfoil  shapes  or  boundary  layer  profiles,  it  is  probably 
fair  to  assume  that,  in  general,  the  energy  in  secondary  flows  for  this  type 
of  airfoil-tunnel  wall  intersection  is 


4 


4 
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D*  =  K1(t/c)4f(n)  .  (4) 

If  the  function  f(n)  in  Fig.  12  is  linearized  over  the  portion  of  the  curve 
* 

0.0  <  6  /c  <  0.1,  then 

f ( n)  =  (f(n)max/(0.1))(6*/c)  =  K^fi/c)  .  (5) 

* 

In  Eq .  (5),  the  displacement  thickness,  6  ,  has  been  assumed  proportional 

* 

to  the  more  frequently  documented  boundary  layer  thickness,  6.  Thus,  D 

e 

becomes 

D*  =  K3(t/c)4(6/c)  .  (6) 

Functionally,  the  drag  correction  was  assumed  to  take  the  following 

form: 


Acd  =  Sl^De’  V  V  ^  *  (7) 

The  inverse  relationship  between  Ac^  and  Dg  has  been  established  by  Barber^®. 
Therefore,  Eq.  (7)  can  be  written  as 


Acd  =  K4 


Mv  V  “» 

(6/c)(t/c)4 


The  effects  of  c„,  c^,  a  and  AR  were  derived  empirically  from  the  available 
experimental  data. 

The  required  Ac^  correction  to  the  drag  data  is  shown  by  the  lower  curve 
in  Fig.  6.  The  deviation  of  drag  is  essentially  zero  at  c  =  0.0.  The  Ac 

d 

curve  increases  to  the  region  where  c^  is  no  longer  constant  and  then 

a 
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decreases.  If  c^  represents  the  slope  of  the  lift  curve  in  the  linear 
ao 

portion  of  the  c^  vs  a  curve,  then  the  shape  of  the  Ac^  vs  a  curve  seems  to 


vary  as  [c^/c^] 


1/2 


It  was  assumed  that  this  slope  variation  represented 


the  angle  of  attack,  a,  contribution  to  the  drag  correction  so  that  Ac  ^  can 
be  expressed  as 


( c£  /c£  )1/2SnK>  cd*  AR) 

a  a 

o 


Acd  =  K5 


’3 ^  l'  d’ 


(6/c) (t/c) 


(9) 


As 


the  experimental  data  show,  Ac^  increases  directly  with  c^;  and  Ac^  is  zero 


at  c^  =  0.0.  This  also  implies  that  the  balance  measures  the  correct  value 

of  c^  at  c^  =  0.0,  namely,  c^  .  Thus,  c^  was  included  as  the  c^  term  which 

o  o 

’’individualizes"  the  correction  procedure  to  specific  airfoils.  When  the 


linear  dependece  on  c  and  the  c,  are  introduced  in  Eq .  (9),  Ac,  becomes 

id  d 

o 

•  cd  K  lci  )‘/2 

o  a  a 

Ac,  =  K  - - S. (AR)  . 

d  o  r  *  i  \  r .  I  i4  4 


(6/c)(t/c) 


(10) 


The  effect  of  aspect  ratio  was  approximated  by 


g4(AR)  =  (AR)' 


(11) 


The  values  of  and  b  were  determined  empirically  from  the  experimental 

data.  It  was  found  that  the  best  fit  to  the  data  was  obtained  for 

K..  =  1.9  x  10  5  and  b  =  -  1/2.  With  the  evaluation  of  these  two  constants 
6 

the  final  form  of  the  equation  to  correct  the  balance  measured  drag  data  for 
the  effect  of  the  sidewall  boundary  layer  is 
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(%)(cd  )(c4  fzi  ) 


1/2 


Ac,  =  1.9  x  10 
a 


-5 


a  a 


(6/c)(t/c)4(AR}1/2 


(12) 


In  summary,  the  proposed  correction  to  balance-measured  drag  data  for 
the  effect  of  the  sidewall  boundary  layer  is  a  function  of  the  airfoil/hydro¬ 
foil  geometry  (t/c  and  AR) ,  the  thickness  of  the  sidewall  boundary  layer 

(6/c),  and  the  accurately  measured  balance  data  fc,  and  c  vs  a). 

v  d  Z  J 

o 


Application  and  Discussion 

Shown  in  Fig.  11  are  the  results  of  applying  the  correction  procedure  to 
the  ARL/PSU  data.  As  can  be  seen  in  this  figure,  the  corrected  drag  polar  is 
in  good  agreement  with  the  ARL/PSU  reference  data. 

In  practice,  this  correction  for  the  sidewall  boundary  layer,  Eq.  (12), 
is  applied  to  the  drag  data  after  the  tare  readings  and  the  traditional 
corrections  have  been  applied.  The  technique  will  be  illustrated  by  applying 
the  procedure  to  data  measured  by  Ward*^  at  CIT  for  the  Canadian  Defense 
Research  Establishment  Atlantic  (DREA).  These  data  were  located  by  the 
authors  after  the  completion  of  Jacobs’  initial  studies  and  were  not  included 
in  the  development  of  Eq .  (12). 

Ward  used  a  three  component,  mechanical  balance  to  measure  lift, 
drag  and  pitching  moment  on  a  6.0  in.  x  6.0  in.  (152.4  mm  x  152.4  mm) 

NACA  16-309  hydrofoil  in  the  CIT  High  Speed  Water  Tunnel.  The  resulting 
noncavltatlng  data  at  50.0  fps  (15.24  m/s)  with  the  tare  corrections 
included  are  shown  by  the  open  circles  in  Figs.  12  and  13.  The  data  were 
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further  corrected  according  to  Pope8  for  solid  blockage,  wake  blockage  and 
lift  effect  (streamline  curvature).  The  sidewall  of  the  test  section  was 
adjusted  to  eliminate  the  horizontal  buoyancy  effect.  The  results  of 
applying  these  traditional  corrections  to  the  data  are  shown  by  the  open 
squares  in  Figs.  12  and  13.  Also  shown  as  solid  lines  in  these  figures 
are  the  NACA  reference  data  (after  correction  for  compressibility  effects) 
measured  by  Lindsey,  et  al20  at  a  Mach  number  of  0.3.  The  required  Acd 
correction  for  Ward's  data  is  shown  by  the  dash-dot  curve  in  Fig.  13. 

Again  it  is  interesting  to  note  that  the  balance  has  measured  the 
correct  drag  value  at  the  zero  lift  condition. 

From  Ref.  [11]  the  following  parameters  were  obtained  for  the 
application  of  Eq .  (12): 

cd  =  0.0009  AR  =  1.00  (t/c)  =  0.09  . 

o 

The  c £  term  is,  of  course,  the  corrected  lift  coefficient  (open  squares)  at 

1/2 

each  a.  The  (c^  / c^  )  term  was  computed  by  fitting  a  differentiable 

a  ao 

mathematical  spline  curve  through  the  corrected  c^  vs.  a  data.  Ward^1 
documents  the  test  section  boundary  layer  characteristics  at  the  balance 
shaft  location  of  the  HSWT.  The  average  value  of  5/c  at  30.0  fps 
(13.24  m/s)  is  0.123.  With  these  values  F,q .  (12)  can  now  be  evaluated. 

The  resulting  Ac  values  are  shown  as  the  solid  squares  in  Fig.  13.  When 

a 

these  computed  Ac^  values  are  applied  to  Ward's  drag  polar,  the  corrected 
data  are  shown  by  the  solid  circles.  The  agreement  is  certainly  encouraging. 
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Conclusions,  Limitations  and  Recommendations 

The  results  of  this  ARL/PSU  experimental  investigation  and  subsequent 
data  analysis  have  revealed  or  reaffirmed  several  important  conclusions 
relative  to  balance-oriented  two-dimensional  airfoil/hydrofoil  testing: 

(1)  The  effect  of  a  small  end  gap  between  the  airfoil  tip  and  the  channel 
wall  is  negligible  provided  the  gap-to-chord  ratio  is  <  0.002. 

(2)  For  gap-to-chord  ratios  <  0.001,  the  effects  of  flow  in  the  region 
of  the  supporting  shaft  are  negligible. 

(3)  Uith  the  application  of  only  traditional  and  tare  corrections,  valid 
drag  polars  can  he  obtained  by  combining  the  balance-measured  values  and 
c^  data  from  wake  traverses. 

(4)  The  disagreement  between  the  traditionally  corrected  balance-measured 
drag  data  and  the  reference  values  is  primarily  the  result  of  the  inter¬ 
action  of  the  airf oil/hydrofoil  and  the  sidewall  boundary  layer. 

(5)  The  effect  of  the  sidewall  boundary  liiyer  on  balance  measured  drag  data 
can  be  accounted  for  by  the  application  of  Eq .  (12). 

The  previous  conclusions  are  not  without  some  limitations.  The  empirical 

development  of  Eq .  (12)  was  conducted  in  the  absence  of  data  obtained  from 

studies  in  which  there  was  a  significant  variation  in  aspect  ratio  or 

thickness  to  chord  ratio.  The  linearized  adaptation  of  Hawthorne's  f(n) 

*  * 

curve  is  only  valid  for  o  /c  <  0.1.  For  values  of  5  /c  >  0.1,  a  different 
approximation  of  f(n)  would  be  required. 
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• 

For  further  refinement  of  Eq .  (12),  one  of  the  recommendations  for  future 
study  would  be  the  acquisition  of  a  larger  data  base,  particularly  with  • 

respect  to  AR  and  t/c  variation.  The  effects  of  the  sidewall  boundary  layer 
should  be  further  investigated  by  wind  tunnel  tests  in  which  the  horseshoe 
vortex  can  be  removed  by  suction  or  blowing.  And,  of  course,  an  attempt  • 

should  be  made  to  extend  the  proposed  correction  procedure  to  hydrofoils 
operating  in  the  cavitating  flow  regimes. 

• 

t 


» 
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Figure  1.  The  two-component  force  balance  which  rotates 
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NACA  0012  AIRFOIL 
BALANCE  DATA 

JACOBS10,  ARL/PSU,  1980 
Re  =  330,000 

□  A?  =  2.04;  C  =  9.0" (22.86mm) 
a  /R  =  1.02;  C  =  18.0" (43. 72mm) 
TRADITIONAL  CORRECTIONS  APPLIED 

—  OHIO  STATE  UNIVERSITY  COMPUTATIONS 
Re  =  330,000 

—  LOFTIN  AND  SMITH1  MEASUREMENTS 
Re  =  700,000 


ANGLE  OF  ATTACK:  a(deg) 


Figure  4.  Sectional  lift  characteristics  measured  by  the 
two-component  force  balance. 


P0  :1N3 10 1 JJ303  ovaa  1VNOI103S 


SECTIONAL  LIFT  COEFFICIENTS 


SECTIONAL  DRAG  COEFFICIENT; C 


0.05 


0.04 


0.03 


0.02 


0.01 


3  March  1984 
ALT : PPJ :GBG : lhz 


0.06  h  NACA  0012  DRAG  POURS 


JACOBS  ,  ARL/PSU 
Re  =  330,000 
REFERENCE  DATA  FIG.  5 

iC i  FROM  BALANCE  MEASUREMENTS 

Cd  FROM  WAKE  SURVEYS 

fR  =  2.04,  C  =  9.0"  (22.86mm) 

BALANCE  MEASUREMENTS  WITH 
TRADITIONAL  CORRECTIONS  (Cg  AND  C.) 

— o-R  =  2.04:  C  =  9.0"  (22.86mm) 
-a-/R  =  1.02;  C  =  18.0"  (45.72mm) 

REQUIRED  DRAG  CORRECTION  V 


"7  A  C 


0.2  0.4 

SECTIONAL  LIFT  COEFFI  Cl  ENT;  C 


Figure  6.  NACA  0012  drag  polar  measured  by  Jacobs  a 
the  required  correction,  Ac,  ,  to  the  drag 
data . 
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WAKE  SURVEYS  AT  VARIOUS  SPAN  LOCATIONS 
NACA  0012  AIRFOIL 

JACOBS10,  ARL/PSU 

a  =  8.4°,  Re  =  300,000,  C  =  9.00'  (22.86mm) 

—  1.75'  (4.45mm)  FROM  BALANCE  END 
••••  1.875'  (4.76mm)  FROM  BALANCE  END 

—  2.0'  (5.08mm)  FROM  BALANCE  END 

—  9.0'  (22.86mm)  FROM  BALANCE  END  (MID-SPAN)  / 


•4  1 _ , 

0 


0.90 


0.95 


DYNAMIC  PRESSURE  RATIO, 


1.00 


Figure  7.  Wake  profiles  in  the  vicinity  of  the  airfoil-wall 
interaction  showing  the  growth  of  a  second  wake 
behind  the  airfoil's  upper  surface. 
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Figure  11.  Jacobs'  (ARL/PSU)  balance  measurements  corrected 
by  Eq .  ( 12) . 


(TARE  CORRECTIONS  ONLY) 


in t a  measured  by  Ward  at  CIT  with  a  6.0  in.  (152.4mm) 
hord  NACA  16-509  hydrofoil  in  the  High  Speed  Water 
untie  I  . 
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